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AN EXPERIMENTAL STUDY OF THE VALIDITY OF THE 


HEAT-FIELD CONCEPT FOR 
SONIC -BOOM ALLEVIATION 

by Rudolph J. Swigart 
The Aerospace Corporation 


SUMMARY 


An experimental program was carried out in the NASA- Langley 4 ft x 
4 ft supersonic pressure tunnel to investigate the validity of the heat-field 
concept for sonic-boom alleviation. The basic idea of the heat-field concept 
is to heat the flow about a supersonic aircraft in such a manner as to obtain 
an increase in effective aircraft length and yield an effective aircraft shape 
(area distribution) that will result in a shock-free pressure signature on the 
ground. By eliminating the shock waves in the classical N-wave sonic-boom 
signature, a substantial reduction in sonic-boom noise heard out-of-doors 
can be achieved. 

In the experimental program, seven l/360-scale wind-tunnel models 
were tested to provide a step-by-step^ verification of the heat-field concept. 
First, a ba^ic body-of-revolution with its lift equivalence in volume 
was tested to provide a baseline pressure signature containing front and rear 
shock waves. Second, a model having- a 5/2-power area distribution was 
tested which, according to theory, should yield,a linear pressure rise with no 
front shockwave. Third, the concept of providing. the 5/2-power area distri- 
bution lay using either an off-axis slender fin below the basic body or an addi- 
tional body of revolution, area-ruled such that the total areas of the basic 
body and the fin or additional body have the. same Mach-plane area distribution 
as the 5/2-power spike, was investigated. Then, to investigate the use of a 
heat field for providing effective body shapes in supersonic flow, a substantial 
portion (approximately 40 percent) of the solid fin was replaced by a heat field 
generated by passing heated nitrogen through the rear of the fin. 

In addition to the above tests, the ability of the off-axis fin to produce a 
finite rise -time signature in the presence of lift was investigated by testing a 
wing-body model both with and without the solid fin previously tested in con- 
junction with the basic body. 

The measured pressure signatures agreed well with theoretical predic- 
tions, with the exception of the signatures for the dual body-of-re volution 
model area-ruled to have the same signature as the 5/2-power frontal spike. 
This model exhibited a strong compression off the nose of the upper body that 
could not be eliminated by compensating expansions from the lower body. 
Comparisons of pressure signatures for the thermal-fin model with no flow, 
cold flow, and heated flow through the fin demonstrated substantial effects due 
to flow and due to heat, essentially verifying the use of mass flow and heat 
fields for providing equivalent area in supersonic flow. Morever, the results 
for the wing -body with solid fin model clearly demonstrated the effectiveness of 
the fin in producing a finite rise -time signature in the presence of lift. 



INTRODUCTION 


An aircraft In supersonic flight generates a complex pressure field that 
simplifies as it propagates, giving rise to the so-called sonic boom (fig. 1). 
The far -field signature is typically N-shaped, being characterized by two 
shockwaves separated by a linear expansion. For cur rent -generation air - 
craft of the supersonic transport (SST ) class, it is estimated that the reflected 
shock strength will be on the order of 100-200 N/m^ (2-4 psf) during cruise 
conditions. ' ‘ 


The principal aspects of the sonic boom causing public concern are 
structural damage and noise. Based on extensive studies of the structural 
damage question, relatively few structures will be damaged by booms having 
overpressures below about 250 N/m^ (5 psf) (ref. 1). Thus, it appears that 
adverse public reaction is associated mainly with the acoustic effects of the 
boom. A substantial startle effect is produced by a shock having an over- 


pressure of the order of 100 N/m^ 
ance. 


or more, and this gives rise to annoy- 


The question, however, of what constitutes an acceptable sonic-boom 
pressure signature is unanswered at this time (ref. 2). It is known that re- 
placement of the front and rear shock waves in the class icaTN- wave sonic- 
boom signature by pressure increments occurring over times as small as 
\ 10 msec results in a substantial reduction of acoustic power in the frequencies 
to which the ear is most sensitive, thereby reducing the outdoor annoyance of 
the boom (refs. 1 and 3). Studies to determine the optimal configuration and 
minimum length for an aircraft having a signature without shock waves have 
been undertaken by"a number of investigators over the past several years 
(refs. 2, 4, 5, 6) with a closed-form result evolving only recently (refs. 2,4). 
These results yield a required length of 173. 74 m (5 70 ft) for an SST -type 
aircraft with a gross weight of 272,155 kg (600,000 Ibm) flying at an altitude of 
18,288 m (60,000 ft) at Mach 2. 7 .- Although this length is substantially below 
earlier estimates of about 304. 8 m (1000 ft) (ref. 6), it is still beyond that 
generally considered feasible for practical implementation. Consequently, 
consideration has been given to the possible use of heat or force fields distri- 
buted about the aircraft in such a manner as to create a so-called "phantom 
body" having the length and area distribution required to produce a bangless - 
boom signature. 

Several investigators have examined heat and force -field concepts and 
have made estimates of the power requirements and weight penalties asso- 
ciated with their implementation (refs. 1, 7, 8, 9). Cheng and Goldburg (ref. 7), 
upon close examination of the use of electro-aerodynamic techniques to re- 
duce the sonic boom as first proposed by Cahn and Andrew (ref. 10), con- 
cluded that on the order of 454,000 kg (1,000,000 lb) of electrical equipment 
would be required to generate the power needed to reduce the overpressure 
of the leaciing shock wave in the sonic-boom signature by ten percent. Both 
Batdorf (ref. 1) and Miller and Carlson (ref. 8), considering a heat field 
generated by direct combustion of the same fuel as used by the main SST en- 
gines, concluded that the front shock wave could be eliminated through properly 
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controlled burning of an amount of fuel corresponding to about 20 percent of 
that normally carried to meet propulsion needs. Oh’ consideration of rear 
shock elimination also, Miller and Carlson concluded that this penalty would 
increase to about 60 percent additional on-board fuel with no account taken of 
the weight and performance penalties of the system. Lipfert (ref. 9) predicted 
the weight penalty associated with alleviation devices in terms of additional 
structural weight and additional fuel weight for both shock-wave elimination 
and to overcome the additional drag due to the additional weight. He considered 
the use of hydrogen (a pollution-free fuel) as well as hydrocarbon fuels not 
only for shock elimination but also to supply the propulsive power to over- 
come the drag penalty. The minimum weight penalty to remove both the front 
and rear shocks for 272,155 kg (600,000 Ibm) aircraft flying a 5000-second 
mission at Mach 2. 7 at an altitude of 18,898 m (62,000 ft) was estimated to be 
226. 796 kg (500,000 Ibm). These tentative conclusions indicate a questionable 
practicality of using free combustion to remove both front and rear shocks. 
However, the use of free combustion in conjunction with other ideas, such as 
the use of a high tail to obtain, in the area-rule sense, the required additional 
length needed to remove the rear shockwave (e.g., ref. 11), could bring such 
concepts back into the realm of practicality. In implementation of these con- 
cepts employment of the area rule allows off-axis combustion to be considered 
(refs. 1,8). This has particular attraction in the case of the front shock, 
since it replaces the difficult problem of burning against the flow ahead of the 
aircraft with one of burning in the direction of flow beneath the aircraft. 

Questions concerned with the engineering and economic feasibility of 
practical incorporation of the previously discussed design features in super- 
sonic transport aircraft can be answered only after exhaustive studies by 
qualified design teams with representatives from all technology areas. How- 
ever, there are rather fundamental questions concerning the validity of the 
heat-field concept than can be answered by means of properly conducted but 
relatively simple wind tunnel tests. Such is the purpose of the present inves- 
tigation. To investigate the feasibility of using an off-axis heat field for re- 
moval of the leading shock wave in the sonic-boom s ignature , an experimental 
program was carried out in the NASA-Langley 4 ft x 4 ft supersonic pressure 
tunnel. Seven wind-tunnel models were tested and pressure signatures re- 
corded at two and three body lengths in the plane directly below the aircraft . 

The results of these tests are reported herein. 

. The author would like to acknowledge the significant contributions of his 
colleagues at The Aerospace Corporation to this program. Dr. S. B. Batdorf 
provided overall supervision to the program and originated the basic concept 
that was investigated, i.e., the use of heat in an off-axis configuration such as 
a thermal fin; Mr. W. R. Grabowsky carried out the thermal analysis and pre- 
liminary design of the heating system incorporated in the thermal fin, and 
Mr. D. A. Durran transformed the preliminary design of all the models into 
detail designs and designed the power -supply console. All of the wind tunnel 
models were fabricated by the Laboratory Operations Machine Shop. 
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SYMBOLS 


Dimensions of physical quantities are given in the International System 
of Units (SI). Equivalent values are indicated parenthetically in the U. S. 
Customary System. Details concerning the use of SI, together with physical 
constants and conversion factors, are given in reference 12. 

A area 


AR 

^L 

c 

P . 
D 

E 

F(y) 


h 

h(y) 

I 

K 

K 

P 

K 

V 

k 

L 


M 


m 


aspect ratio 
lift coefficient 

specific heat at constant pressure 
heating -tube diameter 

elliptic integral of the second kind 

1 / A" (^) 

Whitham F function, — * — ^ 


zlrf 


v/ y 




where A”(|) is the second derivative of the cross-sectional 
area distribution 

static enthalpy; also convective heat-transfer coefficient . 

fin width in maximum -thickne s s plane 

current 

constant in 5/2-power law body -shape equation 

/ 2 

potential lift coefficient/sinOtcos Ct 

2 

vortej; lift coefficient /sin acosce 
?+l 

~ also thermarl conductivity 

longitudinal separation distance of bodies in twin-body 
model (fig, 6) 

various fin Mach -plane chord lengths, depending on subscript 
(Appendix A) 

Mach number, also me /pD 
■ P 

mass -flow rate; also parameter of elliptic integral of 
second kind 
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M. W. 

molecular weight 

P 

2 . 

I R/TTDt, electrical power per unit area of heating-tube wall 

P 

static pressure 

Ap 

difference between disturbance pressure caused by model and 
freestream static pressure 

Pr 

Prandtl number 

'ic 

convective heat -transfer rate 

R 

resistivity 

r 

perpendicular distance from model to measuring probe; also 
heater-tube radial variable 

Re 

Reynolds number 

S 

perpendicular separation distance of bodies in twin -body 
model (fig. 6); also heating-tube cros s -sectional area 

T . 

static temperature 

T° 

stagnation temperature 

t 

heating -tube wall thickness 

V 

velocity 

w - 

fin trailing -edge width (Appendix A); also nozzle width 
(Appendix B) 

X 

distance measured along longitudinal axis from model nose; 
also distance along heating tube 

y 

vertical distance measured from tip of fin; also characteristic 
variable in Whitham theory (Appendix A) 

a 

angle of attack; also resistivity temperature-dependence 
coefficient 

ys 

/ 2 

VM - 1; also coefficient of temperature dependence of 
thermal conductivity 

y 

specific -heat ratio 

€ 

emissivity 

^LE 

leading -edge sweep angle 
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^ Mach angle, sin“l ^ 

^ . M 

P density 

<7 Stefan -Boltzmann constant 

Subs cr ipts 

B basic body of revolution 

b bulk 


c 

D 

e 

F 

H 

j 

L 

o 

R 

S 

w 

P 


centerline value 

based on heating -tube diameter 

effective or exit 

fin structure 

heat 

properties in nozzle jet at exit plane 
due to lift 

initial or entrance value 
reference condition 
5/2 -power frontal spike 
wall condition 
in Mach planes 


GO 


free -stream .value 



1. TEST FACILITY 


The tests were carried out in the NASA-Langley 4 ft x 4 ft supersonic 
pressure tunnel at a free -stream Mach number of 2. 01 and total temperature 
of 317°K (570°R). All of the tests were performed at a tunnel stagnation pres- 
sure of 69 kN/m^ (10 psia) with the exception of two signatures for the ther- 
mal-fin model that were taken at a tunnel stagnation pressure of 51. 7kN/m^ 

(7. 5 psia). Many tests in which sonic-boom pressure signatures have been 
measured for both simple and complex configurations have been carried out 
in this facility, and high-resolution signatures have been obtained (refs. 13, 
14,15). 

For the tests under consideration, the models were mounted on a re - 
motely controlled actuator permitting their movement relative to a slender 
static-pressure probe (fig. 2). To obtain a signature, the probe position is 
held fixed in the tunnel and the model is moved forward relative to it. In this 
manner, the flow disturbances originating at different positions on the body 
always traverse the same region of flow between the body and probe, thereby 
having any free -stream flow nonuniformities in common. In this fashion, dif- 
ferent points in the signature are referenced to a common base. 

2. TEST MODELS AND MODEL DESIGN 


Seven wind-tunnel models were tested during the course of the investi- 
gation. Detailed design rationale and mathematical description of the models 
are given in Appendix A. A basic body-of-revolution (fig. 3) having a nominal 
length of 25.4 cm (10 in. ) was used to represent at 1/360 scale a typical SST 
configuration. Such a body provides a baseline signature containing front and 
rear shock waves. The basic-body area development was patterned on an SST 
area distribution taken from reference 6 and was modified slightly so that the 
front shock wave could be eliminated by heat addition only (i.e., no heat 
removal) (fig. 4). In this model, the lift is replaced by its equivalence in 
volume according to the rules of linearized theory (ref. 16). In addition, the 
nose of the model is modified to have a 5/2-power area distribution such that 
the front shock wave will not form within a distance of about 25 cm (10 in. ) 
from the body (Appendix A). The reason for this is to allow the fin of the 
model discussed below to interact with.the basic-body disturbances before 
nonlinear effects become important, thereby satisfying requirements for ap- 
plicability of the supersonic area rule. 

The second model tested is 'a body of revolution whose area grows as the 
5/2-power distance from its nose (fig. 5). This is the minimum-length area 
distribution for a finite rise-time pressure signature, and is the distribution 
that will result in a linear pressure rise for the first compression in the sig- 
nature according to the basic theory of Whitham (ref. 17). This model was 
designed to represent an airplane equivalent body with a full scale length of 
152 m (500 ft), and be 61 m (200 ft) longer than the airplane itself (see fig. 6), 
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Figure 3. Basic Body of Revolution Representing 
Typical SST Configuration 










and of such a contour that the area difference between the basic body and the 
5/2-power frontal spike is constant aft of the nose of the basic body. This 
would correspond to the case of heat distributed upstream of the basic body in 
such a fashion as to displace streamlines in a manner identical to the 5/2- 
power frontal spike. The heating would end at the nose of the basic body 
(ref. 1). 

The third model incorporates a solid fin situated below the basic body of 
revolution and is designed so that the area distribution in Mach cutting planes 
perpendicular to the lateral symmetry plane of the fin is identical to that of the 
5/2-power spike (fig. 7). The purpose of the 5/2-power area distribution nose 
on the basic body (discussed above) is to allow disturbances from the fin to 
interact with those from the basic body before nonlinear effects become im- 
portant, thereby satisfying requirements for applicability of the supersonic 
area rule. The fin is designed to have a double -wedge cross-section in Mach 
planes, with a maximum front-wedge angle of four degrees in horizontal 
planes. It is necessary to keep the front wedge-angle in horizontal planes 
small in order to remain within the bounds of linearized theory, upon which 
the supersonic area rule is based (ref. 16). When these conditions are met, 
the pressure signature below the fin in the lateral symmetry plane of the fin 
is theoretically identical to that of the 5/2-power spike. 

The fourth model tested is the so-called twin -body wherein a 
second body of revolution is situated below the basic body of revolution (fig. 8). 
The second body has a Mach-plane area distribution defined so that the sum of 
its area and that of the basic body is equal to the area distribution of the 5/2- 
power spike model. Again, according to the area rule, the pressure signa- 
tures of this model in the plane through both axes of the bodies of revolution 
should be identical to that of the 5/2-power body. The model is designed for 
separation distances, S, between the two bodies of 4. 445 cm (1.75 in. ) and 
8. 89 cm (3. 50 in. ) and variable relative length distances, L. 

The fifth model tested is the thermal-fin model in which approximately 
40 percent of the fin volume of Model 3 is replaced by a heat field designed to 
have the equivalent area distribution in Mach planes of the solid fin (fig. 9). 
The heat field is generated by passing nitrogen through a 0.3175-cm (1/8-in. ) 
diameter heating tube that continues through the model into a nozzle whose 
exit plane is at the rear of the fin (fig. 10). In order to accommodate the 
1/8-in. heating tube in the fin section, the scale of the thermal-fin model was 
increased by a factor of 1.5 over its solid counterpart. Hence, the basic- 
body portion of the thermal fin model has a nominal length of 3 8. 1 cm (15 in. ). 

The thermal fin Mach-plane area distribution is shown in figure 11. In 
the Mach plane at the end of the heated region, better than 60 percent'of the 
fin area is represented by heat. This percentage drops to 43 percent at 
30. 163 cm (11. 875 in. ) and remains constant thereafter. In terms of volume, 
about 40 percent of the total fin volume is represented by the heat field. 

Details of the design of the thermal fin and heating tube are given in 
Appendices B and C. 
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Figure 9. Thermal Fin Model 





NICHROME TUBE 



Figure 10. Thermal-Fin Model Cross Section 





Figure 11. Thermal-Fin Mach-Plane Area Distribution 





The sixth model tested is the so-called wing-body model consisting of a 
body of revolution and a double delta wing for which the sum of the Mach 
cutting-plane areas due-to-volume and the equivalent Mach-plane area due -to- 
wing lift is equal to the Mach-plane area distribution of the basic body of revo- 
lution (fig. 12). Then, according to linearized theory and the basic sonic- 
boom signature theory of Whitham, this model should have a pressure signa- 
ture below it in the lateral symmetry plane through the fuselage centerline 
equal to that of the basic body. 

The wing, is a 67. 5-degree sweep-angle double delta set at 40 minutes 
angle of attack. The lift for the wing was calculated using the theory of 
Polhamus (ref. 18); details are given in Appendix A. At the end of the body, 
the equivalent area due'-to-lift represents about 30 percent of the total effec- 
tive area. 

The seventh model tested is the wing-body model with the solid fin of 
Model 3 attached below it. The purpose of this model was to test the effective- 
ness of the fin in generating a finite rise -time signature in the presence of a 
lateral lift distribution.' 


3. RESULTS AND DISCUSSION 


The results for the measured pressure signatures are shown in figures 
13-29. In these figures, the nondimens ionalized overpressure 4p/poo >-s 
shown as a fxmction of distance, x - fix , behind the Mach line through the 
point of initial model disturbance at a fixed perpendicular distance r frorn the 
centerline of the body (fig. 2). Signatures were taken at values of r = 

50. 8 cm and 76. 2 cm (20 in. and 30 in. ) for all except the thermal-fin model, 
corresponding, respectively, to two and three body lengths for the nominal 
ten-inch models. For the 38. 1-cm (15-in. ) thermal-fin model, signatures 
were taken at r = 76.2 cm (30 in. ), corresponding to two body lengths. Thus, 
all signatures are of the near field. 

Figures 13 and 14 show the results for the basic body-of-re volution 
airplane. Cornparison is shown with theoretical predictions made using a 
computer code developed at NASA-Langley (ref. 13) based on the theory of 
Whitham (ref. 17). Note that, in general, the agreement with theory is good. 
The magnitude of the front shock overpressure and rear expansion are under - 
predicted, however, and the positions of the rear expansion and rear shock 
are somewhat in error, A possible explanation for the latter is the fact that 
the theory assumes that all disturbances originate oh the body axis rather than 
on the body surface, which would cause the theoretical disturbance to be dis- 
placed downstream. 

Figures 15 and 16 show a comparison of theoretical and experimental 
pressure signatures for the frontal spike model. Note the excellent agree- 
ment with theory in this case, the measured data very nearly following the 
predicted linear pressure rise. 
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Pressure Distributions 




FRONTAL SPIKE, r = 50.8 cm (20 in.] 



Pressure Distributions 






FRONTAL SPIKE, r = 76.2 cm (30 in.) 
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The most important data, of course, are those generated by the solid 
fin, the twin-body, and the thermal-fin models. The results for the solid fin 
are shown in figures 17 and 18. Recall that the fin is area-ruled such that 
the area distribution in Mach cutting planes is identical to that of the 5/2- 
power frontal spike. In addition, the fin -is designed such that disturbances 
to the free stream are small. ■ When these two conditions are met, a pressure 
signature in the plane of symmetry below the fin identical to that of the 5/2- 
power spike is predicted. Note that such a result is nearly obtained in 
figures 17 and 18. The important question, of course, is whether any local 
compression regions will cause a shock wave to develop at larger distances. 
Two such regions are seen at r = 50. 8 cm (20 in. ) in figure 17, one at the 
beginning and one' toward the end of the linear pressure-rise region. When ^ 
the signature at r = 76.2 cm (30 in. ) in figure 18 is considered, however, it 
is seen that no significant steepening has occurred for the compressive region 
at the beginning, and the one toward the end of the linear region has nearly 
disappeared. Even if shockwaves developed from these compressions, how- 
ever, they would be no more than ten percent of the total pressure rise. 

The results for the twin-body model were not as encouraging, however. 
These are shown in figures 19, 20, 21 and 22. Figures 19, 20 and 21 give the 
signatures at r = 50. 8 cm (20 in. ) for a separation distance, S, of 8. 89 cm 
(3.5 in. ) arid various settings of relative lengths, L (see fig. 8 for definitions 
of S and L). The Lq corresponds to the nominal setting as would be predicted 
by linearized theory, i.e., all Mach lines at their undisturbed free-stream 
positions. Figure 19 shows a local steep compression occurring about 15 cm 
into the signature, followed by an overexpansion and subsequent re compres- 
sion. This is caused by the disturbance coming off of the nose of the basic- - 
body of the twin-body model, the signature up to that point being due only to 
the lower body. Recall that the basic body was designed with a local 5/2- 
power cusped nose 1. 359 cm (0. 535 in. ) in length such that a shock would not 
form prior, to a radial distance, r, of about 25 cm (Appendix A). Expansions 
from the shoulder of the lower body are designed to interact and partially 
cancel the compressions originating from the cusped nose, thereby maintain- 
ing the linear growth. Since the Mach number of the tunnel flow is in actuality 
slightly above two (2. 01), the first relative -length setting, L, was made at a 
value 0. 1588 cm (0. 0625 in. ) less than nominal to account for this. An initial 
compression was noted (fig. 19) and it was thought that perhaps the expansions 
from the lower body were not coming into play soon enough. To correct for 
this, a second run was made with a value of L = 0. 635 cm (0. 25 in. ) ahead of 
nominal (fig. 20). , Note that the initial compression is now replaced by an 
expansion, with a subsequent over compre ss ion and re -expansion. Thus, the 
lower body has been moved too far forward to provide the desired cancellation. 
The results, of moving it back to a position 0. 1588 cm (0. 0625 in. ) ahead of 
nominal are shown in figure 21. Even though the initial expansion is nearly 
eliminated, the strong overcompression is still present. As a matter of fact, 
comparison of figures 19, 20 and 21 reveals that the peak of the local over- 
pressure region remains at the same position and is of the same magnitude 
for all three relative body positions, thereby indicating a negligible effect of 
the lower -body expansion on the compression off of the nose of the upper body. 
One possible explanation for this is the boundary-layer displacement thickness 
effect, which tends to enhance the compression and weaken the expansion. 
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TWIN-BODY MODEL, S = 8.89 ci 
r = 5Q. 8 cm (20 in.), L = Lq - 0. 
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WIN-BODY MODEL, S = 8.89 cm (3.5 in.) 
= 50. 8 cm (20 In. ), L = Lq + 0. 635 cm 
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Comparison of Experimental and Theoretical 
Pressure Distributions 
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This effect could be Investigated by compensating for tKe boundary-layer dis- 
placement thickness in the body contours. The shock wave occurring at about 
X - jffr = 25 cm is due to the compression generated off of the lower body sup- 
port (fig. 8). 

Similar results to those at the 50. 8 cm (20 in. ) station for the 8. 89 cm 
(3-1/2 in. ) separation were obtained at the 76. 2 cm (30 in. ) station. In addi- 
tion, a result for the 4.445 cm (1. 75 in. ) separation is shown in figure 22 at 
the nominal relative -length setting. Results similar to those at the 8. 89 cm 
(3. 5 in. ) separation are noted. 

Thus, based on the more favorable results for the solid-fin model than 
for the twin body, a thermal-fin counterpart was constructed and tested. 

The results are shown in figures 23 and 24 for a wind-tunnel total pres- 
sure of 51. 7 kN/m*^ (7. 5 psia) and in figures 25, 26 and 27 for a total pressure 
of 69 kN/m^ (10 psia). The former is the design-condition setting to obtain the 
proper ratio of jet -to-freestream stagnation pressures as determined by the exit 
static-pressure match condition (Appendix B). Figure 23 shows the resulting 
signature at r = 76. 2 cm (30 in. ) for the case of no flow through the fin nozzle. 
Note the large expansion and subsequent recompression. This portion of the 
signature corresponds to disturbances generated by the nozzle region of the 
fin, the expansion being generated by the growing area deficit through the 
nozzle region, and the recompression caused by this deficit becoming constant 
at the end of that region. Figure 24 shows the result of passing heated nitro- 
gen through the fin at the design mass-flow ra!te and temperature (Appendices 
B and C). Note the substantial reduction in the nozzle deficit region, although 
an expansion still exists at the outset. Visual examination of the model after 
the test indicated a greater nozzle width in the center than at the ends when 
the flow was turned on, giving a mass (and"hence area) deficit in the lower 
region of the fin. This is the probable cause of the expansion-recompression 
region noted in the signature. Further, note the two points, marked ”cold 
flow” and "no flow. " These were obtained by first turning the heater power 
off and then turning the nozzle flow off. In this manner, the effects of the flow 
and heat are separated, and it is seen that they both contribute about equally to 
the overpressure at this point. 

Figures 25,. 26 ahy 27 show the effects of cold and hot flow on the ther- 
mal fin pressure signature for a wind-tunnel total pressure of 69 kN/m^ 

(10 psia). Although the nozzle is underexpanded at this pressure and the effec- 
tive-area contribution is only about three-quarters that of the fully -expanded 
case, note that the effects on the pressure signature are still substantial. 

The final models tested were the wing-body and wing -body with solid fin. 
The objective of these tests was to assess the effectiveness of the solid fin in 
producing a finite rise -time signature when lift is present in actuality as 
opposed to its equivalent area form as in the basic body of revolution. Fig- 
ure 28 shows the results for the wing-body alone compared with the theore- 
tical prediction for the basic body of revolution. Recall that the wing-body 
was designed so that the sum of the areas due to volume of the fuselage and 
wing plus the equivalent area due-to-lift of the wing is equal to the area of the 
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basic body of revolution in each Mach plane. Hence, the wing-body should 
have the same pressure signature as the basic body of revolution. The ex- 
perimental results of figure 27 compare well with the theoretical prediction 
for the basic body and with the experimental results of figure 11, thereby 
' verifying the equivalence. 

Figure 29 shows the results of adding the solid fin previously tested in 
conjunction with. the basic body of revolution (fig. 5) to the wing-body model. 
The resulting pressure signature compares well with the theoretical predic- 
tion except in a small region where the nose of the basic body begins to affect 
the signature. A sudden expansion followed by a gradual re compress ion is 
noted in this region. One probable cause of this is a small gap between the 
fin and basic body in that region that was noted upon detailed examination of 
the model after the test data were taken. 


4. CONCLUSIONS 


Based on the results of the wind-tunnel tests reported herein, several 
conclusions may be drawn. Most importantly, the use of heat fields to pro- 
vide effective areas through streamtube expansion has been verified. In addi- 
tion, it has been shown that off-axis structures and/or heat fields may be used 
to produce finite rise -time sonic-boom signatures. Moreover, such struc- 
tures are effective in producing a finite rise-time signature in the ground 
track in the presence of a lateral lift distribution. 

\ 

5. SUGGESTED FUTURE RESEARCH 


Major criticism of the heat-field approach to sonic-boom alleviation has 
centered upon two areas; first, the lack of this approach to significantly re- 
duce acoustic energy in those frequencies to which structures respond, and 
second the large penalties associated with means of implementing this ap- 
proach into a practical aircraft design. 

The indoor noise characteristics are mainly a function of signature im- 
pulse and overpressure. Since the area distribution for minimum overpres- 
sure varies with distance from the aircraft nose to the 3/2 power whereas 
that for finite rise time has a 5/2 -power variation, it is not possible to 
achieve a minimum overpressure finite rise-time signature. Hence, the best 
that can be done in answer to the first criticism above is to minimize the peak 
overpressure increase associated'with implementation penalties to obtain a 
finite rise-time signature. Cursory studies of these penalties have been 
undertaken by several investigators (refs. 1, 8, 9). A design study of an air- 
craft taking into account the capability of producing a finite rise -time signa- 
ture as an integral part of the design has not been performed, however. Such 
a study would be in order to properly assess and minimize the penalties 
associated with achievement of a finite rise -time signature. 
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WING-BODY WITH SOLID FIN, r = 50.8 cm (20 In.) 
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APPENDIX A 


DESIGN EQUATIONS FOR BASIC BODY, FRONTAL SPIKE, 
SOLID FIN, thermal FIN, AND WING-BODY MODELS 

A. 1 Basic Body of Revolution and Frontal Spike 


The basic body of revolution and 5/2-power frontal spike models were 
designed in conjunction with one another to satisfy the requirement that if the 
5/2 -power frontal spike were simulated by a heat field, the front shock would . 
be eliminated by heat addition only, with the heating ending at the nose of the 
basic body. In addition, in order to obtain a finite rise -time signature on the 
ground in the equivalent full-scale sense for a 272,155-kg (600, 000-lb) air- 
craft flying at 18,288-m (60,000-ft) altitude, it was desired that the frontal 
spike begin 60,96 m (200 ft) in front of the basic body (fig. 6). Considering the 
representative SST area distribution of fig. 4, the 5/2-power contour is main- 
tained up until X - 68.58 m (225 ft). Thus, two equations for the 5/2-power 
constant of proportionality and heated area increment are provided. The re- 
maining contour from 68, 58 m (225 ft) to 91. 44 m (300 ft) was fit with a cubic 
equation. There results, for the model frontal spike area distribution to a 
scale of 1 cm - 3.6 m (1 in. = 30 ft): 


-16. 93 cm < X $ 19. 05 cm (-20/3 in. < x <15/2 in. ) 

Ag = 0. 946 X 10"^ (X + 16. 93)^^^ cm^ 

19. 05 cm < X < 25. 4 cm (15/2 in. <x <10 in.): 

A = -0. 567 X 10’^(x + 16.93)^ + 0.591(x + 16. 93)^ 
o 

- 19. 975(x + 16.93) + 225.47 cm^ 


A 5/2-power cusp was fitted to the nose of the basic body to prevent a 
shock from forming in front of the fin when used with the fin attached (see 
below). To prevent shock formation before a nominal distance of 25.4 cm 
(10 in. ) below the model, it is required that the advance be such that x -/Sr = 0 
at r = 25.4 cm, i.e,, from the theory of Whitham (reL 17), the characteristics 
are given by 

y = X - /Sr + kF(y)r^^^ 



where 


k 


(Y+ 1) 

'/I 



CO 


For a 5/2-power contour, F(y) = 15/^ Ky, where K is the proportionality 
constant to be determine^ i.e., A = is the area growth in terms of the 

axial variable x. Thus, K is given by ^ 


K = 


16 


15 kr 


1/2 


and, for r = 25.4 cm (10 in. J, K = 0.0176 cm - (0. 028 in. ). A value of 

0.y0188 cm“^'^ (0. 03 in. "^'^)^was selected. The length of the cusp is deter- 
mined by matching area and slope on the body to which the cusp is fitted if the 
slope downstream of the matching point is less than that for the 5/2-power- 
contour. For the basic body under consideration, the matching point was 
selected at a scale length of 1. 588 cm (0. 625 in. ), the slope downstream of this 
point being less than for the 5/2-power body. The area of the basic body at 
this point required a total cusp length of 2. 945 cm (1. 16 in. ). Hence, the 
origin of the basic body modified with the 5/2-power nose cusp is at 
X = -1.36 cm (-0. 535 in. ). Thus, the basic body contour is given by 


-1. 359 cm < X < 1. 588 cm (-0. 535 in. < x < 0. 625 in. ): 

A„ = l.:88 X 10"^( X ,+ 1. 359)^^^^ cm^ 

Xj 

1. 588 cm < X £ 3. 175 cm (0. 625 in. £ x £ 1. 25 in. ): 

A^ = 1. 88 X 10'^ (x + 1/359)^'^^ - 1.441 x lO'^ (x - 1. 588)^^^ 

JD 

8. 94 X 10"^ (x - 1. 588)^"^^ cm^ 
3. 175 cm £ X £ 19. 05 cm (1. 25 in. £ x £ 7. 5 in.): 

A^ = 0. 946 X 10"^(x -f 16. 93)^^'^ - 1. 1161 cm^ 

Jj 
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19. 05 cm < X ^ 25. 4 cm (7, 5 in. < x < 10 in. ); 

A„ = -0. 567 X 10“^ (x + 16. 93)^ + 0. 591(x + 16. 93)^ 
Jd ' 

- 19. 975 (x + ,16. 93) + 224. 354 cm^ 


A. 2 Solid Fin 


It is required that the Mach -plane area distribution of the fin plus the 
basic body to which it is attached be equal to that of the 5/2-power frontal 
spike. A double-wedge fin cross-sectional area was selected with straight 
leading and trailing edges for ease of machining. In addition, in order that 
the restrictions of slender -body theory not be violated, a maximum lateral 
flow angle of four degrees was selected. It is also desirable not to allow the ^ 

fin leading edge to be ahead of the body nose in vertical planes at any point. 

Moreover, since the maxiinum cross-sectional width^is growing as distance y 
to the 3/2 power (Sketch A), the maximum lateral flow angle will occur at the 
largest value of y for which straight leading- and trailing -edges are main- 
tained. These considerations require that the straight leading- and trailing - 
edge conditions be abandoned at some point up the fin and the fin width not 
grow thereafter. The point was selected as y = 5. 69 cm (2.25 in. ). All of 
these considerations resulted in the following fin design equations: 

0 ^ y ^ -5. 69 cm (0 < y < 2.25 in. ) : ( 
h(y) = 0. 0195 y^^^ cm 
ij(y) = 0. 476 y • 

= 0. 293 y 

5. 69 cm < y < yj^ (2. 25 in. < y < yj^): 
h(y ) •= 0. 267 cm 

0.4755 y if y < 8. 6l ci 

/i(y) = 

4. 1062 cm if y > 8. 61 ci 


x z 



Sketch A 
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A(y) = 




- /j, y < 6. 769 cm 


/'vA 


2(y- 5. 027). 6, 769 cm<y<y. 


0, y < 6. 769 cm 

4(Ag - Ag) - +^ 2 ) 



, 6. 769 cm < y < y. 


V Sketch B 


where w is the fin trailing-edge width (Sketch B). 


y. < y < 7(^3 = 0) (Sketch C) 


J = 2(y - 5. 027) cm 


f/2 


V 


where h = 0. 267 cm, a = 4. 106 cm and 

^_ and A are obtained as follows. 

^3 seg 

To determine 


Sketch C 



Sketch D 

Mach-Plane Gr os s -Sectional Geometry 



From sketch D 


cos 6 


T^/airvH - - i^) 

_ 


s in ^ = 


I 

R 


therefore 



+ 



/sin// - ily 



R 


2 


( 1 ) 



N 


X Z 




Plane normal to axis 
at X = z - {2^ - l^)cos^ 


Sketch E2 


From sketches El and E2 



( 2 ) 

(3) 


Combining equations (1), (2) and (3) 



2 

r 


■uZ 2 

- ^ cos ^ 

. 2 
sin // 


(4) 
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Also 


taxiQ 


a - i. 


(5) 


Combining equations (1), (4) and (5) 

2 

r (x ) 


l^a - tan^0 - r^(x) + [ ' (>^3 ” = 0 

where a = 4.106 cm (1.617 in.), i. = tan^ = 0.0324 

" 4 cos)U 


^3(2) 


seg 


seg 


2 A / \2 ^ 2 


seg cos^ 



z - (^3 - i^)coSjt/ 

Simpson's Rule: 
b 

y' y dx = |- + 4 (y 2 +74 + -. - + 7 ^^) + 2(73 +75 + 







where b - a = {l^ - i!^)cos//; let n = 4, then h = (i^ - ^^)cos^/ A. 
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Therefore 


I 


(^3 ' -^ 4 ^ r 1 

^seg = 6 + ^4^ +^^3 +^ 5 ] 

Yj =yi^{x.^) - (Xj - x^)^ tan^fl x^ = z - - l^)cosfl 

^2 ^ +h) .= y^z - (^3 - i^)cosfJ j^l - |-jj 


' Therefore 

Y2 ~ (^2 ■ ' ^2 ^ ^ ~ ^^3 ~ (^ " I") 

V C, 2 2 / ■ 1 \ 

r (X 3 ) - (x^ - x^) tan ju x^ = z - (^3 - i^)cos^ Yl - jj 

et cetera. It is now convenient to work in the variable z rather than y. 

The following relationships determine the fin design in the range 
= 0) < z < 25.4 cm (10 in. ):[w(x), determined]. To 

determine 

X = z - r, cotyU (sketch F) 

o b ' ' 


The intersection of the Mach planes and axes are given by 


+ x/3b 


X. - X 

I o 


c. = . y. = 

^ 1 -V^ ^ /T 


where y = ax + b is the equation of Axis I or II, as appropriate. 

Axis I: y = 0. 1763x - 0. 0848 cm 
Axis II: y = 0. 0782x + 0. 214 cm 
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Then 


I2 = ±a/(x.-z)^ + [y- -r^{z)f, .+: X. > z 


X. < z 
L 


a.nd j.yis determined by 

r^(x) - [r|^(z) - - i^)sinfif‘ - 0. 00696 = 0 

where x = z - ^3 - ^2 )cos//. The / is given by 1 = 2(6.215 - y.) cm. Then w 
is determined from 

4(A - AB)-2(hi3-A ) 

w= ^ S-.h 

where h = 0.267 cm and A is determined as before. 

seg 

x z 



c = 2. 54 cm e = 3. 045 cm 

d = 0.362 cm f = 0.4524 cm 


Sketch F. Design Nomenclature 








A. 3 Thermal Fin 


In each Mach plane, the sum of the fin area plus the effective area due- 
to-heat plus the basic-body area must equal the 5/2-power frontal spike area 
distribution, i.e. ' 


where 




+ A 


H 


+ A^ 





fin Mach-plane area due -to-structure 
fin Mach-plane area due-to-heat , 
Mach-plane area due to basic body 
Mach-plane area of 5/2 -power frontal spike 


To accommodate a 0. 3175-cm (l/8-in. ) diameter nichrome heating tube 
inside the thermal fin, it was necessary to make this model one and one -half 
times as large as the solid-fin model. Hence, the basic body portion of the 
thermal-fin model is given by 


-2. 0384 < X < 2. 3813 cm (-0. 8025 < x < 0. 9375 in. ); 

A„(x) = 0. 0154 (x + 2. 0384)^^^ cm^ 

2. 3813 < X < 4. 7625 cm (0. 9375 < x < 1. 875 in. ): 

A^(x) = 0. 0154(x + 2. 0384)^^^ - 0. 0118(x - 2. 3813)^^^ 

- 0. 1091(x - 2. 3813)^^^ cm^ 

4. 7625 < X < 28. 6 cm (1. 875 < x < 11. 25 in. ); 

5/2 2 

A^(x) = 0. 00077(x + 25. 4) ' - 0.9896 cm 

li 

2. 858 < X < 38. 1 cm (11. 25 < x < 15 in. ): 

Ag(x) = -0. 00378{x + 25. 4)^ + 0. 591(x + 10)^ - 29. 962(x + 10) + 504. 799 cm^ 
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As In the case of the solid fin, it is convenient to work in'distance 
variable y along the fin rather than x (see sketch G). These are related by 
y = x/yr. The design equations for the fin are 


0 < y < 8. 5725 cm (0 < y < 3. 375 in. ),: 
h(y) = 0. 0159 y^^^ cm 
i^(y) . 0.4755 y 
/^(y) = 0.2928 y 

The fin heat starts at y = 8. 58 cm (3. 375 
bution to the fin area is 



J 

in. ). From Appendix B this contri- 


A„ = 70.452 w" (y 8. 573) cm^ 

n 

i 

where w''^ is the nozzle throat width. Then, in the region 8. 5 73 < y < 9. 254 cm 
(3. 375 in. < y < 3. 643 in. ): 

^1 

i 

h 
d 

^2 

Then, letting d + d = D there is obtained 

5 /2 

„ Kv ' - 0. 8129 - 0. 1897(y - 8. 573) - (0. 1994 + 2. 73w*)(2. 5103 - 34. 3634w*) - 10. 92w*(v - 8. 573) 

17. 1817-w* - (y - 3. 375) 



where K = 6. 11x10 ^ cm 
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This is of the form D = N/D; N must be 0 when D = 0. This is so for 

w* = 0. 0396 cm 
y = 9. 253 cm 


Then, = 1* 1473 cm. 

The total length of the heated region is given by 


A'!' 


m/p '!-V* 

W^' 


5. 0114 cm 


at the design condition. In 9. 253 cm < y ^ 12, 954 cm (3. 643 in. < y ^ 5. 100 in.) 


^2 = 1* 1473 cm 
II +£3 = 0.4755 y 
h = 0. 399 . cm 
The is given by 

+ ( I + 2. 73w*)i2 

. +10. 92w*(y - 8. 573) = Ky^^^ 


( 



At y - 12.954 cm, 4^3 becomes constant, and at y = 13.487 cm, the basic 
body of revolution begins to contribute to the effective area distribution. 
Hence, in 12. 954 cm < y ^ 13. 487 cm: + £3 - 6. 1595 cm, £2 ~ 1. 148 cm 

and £j^ is determined by the previous expression. 

Then, in 12, 954 cm < y < 13. 584 cm (end of heated region), is deter- 
mined by 

2 A ^ ■'■(2 +2.73w+^i2 ^ ^ 

where h = 0. 399 cm and A is given as before. 
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For y > 13. 584 cm, the heating contribution to the effective area is 
constant at 2, 1697 cm^ (0. 3363 in.^). The I 2 allowed to grow linearly and the 
governing equations are, for 13. 584 cm < y < 13. 754 cm 

j^l + ^3 = 6. 1595 cm 

= 1. 1473 + 2(y - 13. 5834) cm 

d/2 = (2. 5103 - 2017) cm 


\ h{6. 160 - i^) + {j + |)i 2 + .2* 1697 + Ag = Ky 

where h = 0. 399 cm. 


5/2 



At y = 13. 754 cm the leading edge of the fin intersects the basic body, 
and the contribution of the front portion begins to decrease. Thus, in 
13. 754 cm < y< ~ 9) (sketch K) 

' = 3 . 840 cm 

^2 = l*l^'^3+2(y- 13.584)cm 

= ,2. 319 cm 

h =0. 399 cm 


tan^ 




= 0. 0519 


(I + 1)4 + "4 +1^1 - 'A - - A 



where j ^3 and Ag^g are determined as in the case of the solid fin (see pages 
45 through 48. 

' When the point at which ^3 = 0 is reached, it is convenient to work in 

axial variable z rather than normal coordinate y (see page 46), The Mach- 
plane area distribution in the range z (.^3 = 0 )< z < z(l^ = 0 ) is given by 

^2 = 1* 1^”^^ + 2(y - 13. 5834) cm 

= (3.81 - z)/cos/i = (2/^^3)(3.81-z) cm 

h =0. 399 cm 
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The 


Is given by 

2 

^ + [r^(z) - (^5 - l^ysinflf' - r^(x) = 0 

where x = z - (Ir- - If )cosA/. Then w is 

. , DO 

given by 

+ +h4-A + 2. 1697 +A„ = A_ 

\ 2 2 / 2 5 seg 

In z{tc = 0) < z < 15, the fin is designed 
according to the constraints 

w = const = 0.465 cm [its value at z(/^ = 0)] 

= 2[3. 856 - r^(z)] (bottom edge horizontal) 

H - ^ determined by 



^ - r^(x) 


= 0 


where x = z - U-l 2)cos/i 


(h + w) 


i-i 

3. 856 - r, (z) + — = — : 
b ^ 


- A 


seg 


+ 2. 1697 + A„ = A_ 



w 


A, 4 Wing-Body 


It-is required that the wing-body combination with wing at angle of 
attack, Ot, have the same area distribution as the SST representative equiva- 
lent body (fig. 3). 
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Based on the theory of Polhamus (ref. 17), the total lift of a delta wing 
is given by the sum of the potential lift and the vortex lift. Accordingly, the 
lift coefficient is given by 


where 


= K sin 

P 


2 

cos Ot 


K sin ct cosct 

V 


K 

P 


ttAR 

2E 


K 

V 


7T\/[16 - (ARy3)^][AR^ 


16E 


+ 16] 


The AR is the aspect ratio and E the elliptic integral of the second kind of 
parameter m = 1 - (^cotXLj])^, where ^ 1 , and is the leading- 

edge sweep angle. 

A wing having a root chord of l6. 94 cm (6-2/3 in. ) was selected (latter 
two-thirds of body length). It was then necessary to select a sweep angle and 
angle of attack. For this test, a subsonic leading edge was desired to avoid 
attached shock waves from the wing. Hence, X j > 60 deg is required for the 
Mach 2 tunnel flow. The final consideration is that there remain sufficient 
radius of the basic body at the end of the model to allow support by the sting. 
On this basis; a 67. 5 -deg sweep angle and 2/3 -deg angle of attack were 
selected. 


The equivalent area due-to-lift is then calculated as a function of running 
distance, x, from the point of intersection of the wing with the body (sketch M). 
It is given by ' 


(x) 


^-^L(x) = |c^S(x) 


P V 
oo oo 


where S(x) - x tan(7T/2 - X ) = 0.414x for X, = 67. 5 deg. Then, for 
2/3 deg, /3 = /3 LE , LE 


A (x) = 0. 008 X 
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Sketch N 


Then, the design condition determining the body shape is 

■^BV -^WV ^Wh ^ 

where 

A = area of body of wing-body combination due to volume 

= area of wing due -to -volume 

A„,. = area of wing due -to -lift 

W ij 

A^ = basic body of revolution area due to volume 
BB 

2 

Then, defining Agy = TTr and noting other definitions from sketches N and O 
there is obtained 



Sketch O Sketch P 
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r = 


(b + r )tanO 
o 


sin(<pJZ) + tan0cos(f/2) 


. _ 2(b + r )rsin(y/2)- r <f _ 

r + ° -T T,, + 0. 008 

cos 40* BB 


provide two equations for the unknowns r and (f) at each value of x where 
b = X tan 22. 5 deg 

r^ = value of r at beginning of wing - 0. 785 cm (0. 309 In. ) 

tan0 selected at 0. 029 
X - 10/3 

X = 

cos 2/ 3° 
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APPENDIX B 


THERMAL FIN MODEL NOZZLE DESIGN ANALYSIS 


The thermal-fin design in Mach planes aft of the point of initiation of the 
hot gas flow depends on the nozzle design. The nozzle was designed on the 
basis of one -dimensional flow considerations for flow with mass addition. 
Consider flow in a streamtube in planes normal to the fin. 



P + dp 
p + dp 
V + dV 
A + dA 
T + dT 
M + dM 
m + dm 

M. W, + d(M. W. ) 


The conservation equations for the flow through the streamtube are 


m. 


Continuity: 


m 


dp dV dA 
P V A 


(B-1) 


where m. P.V.A.; m = PVA 
J J J J 


Momentunj: (p. - p)A. - Adp mdV + m.(V - V.) 

/ \ 

Energy; m(dh + VdV) + m^ I h - h^ + ^ ^ j ~ ^ 

Stata: ^ 

p P T M. W. 


(B-2) 

(B-3) 

(B-4) 


Combining these equations, yields 


dA 




= :^LJ. 1 + 

A T r I -r 


. y-1 




fM. 


A. 

J 


PV 


M 






P V / ' 

.JU. I _ M. w. 

PV I M. W.. 

\ 


1 + (y-l)M^ /P.i _ A ^ 1-M^ ^ ^ 

7M^ Y / P 


(B-5) 


59 



'y 'y - - “y 

where the relationships V /h = (y-l)M , h°/h =1 +(y-l) Mj /2, h°/h = 

1 +(y-l) M^/2 and d(M. W. )/M. W. = m./m ^(1 - M. W. /M. W..) have been uti- 
lized. ^ ^ 


To avoid shockwaves in the pressure signature, a nozzle design with 
Pj/p - 1, dp = 0 was selected. This reduces equation (B-5) to 


dA 

A. 

J 



- [1 + (y-l)M^] 




M. W. 

M. W. . 
J 


(B-6) 


The P.V./PV can be expressed as 
J J 


P.V. M. 

.. J . J = _1 
P V M 


M. W. 


1 + 


y-1 


M. 


M. W. 


1 


Thus, from equations (B-6) and (B-7) 



dA 

A. 

J 


f 1 , M, 


M. W. . 
1 

M. W. 



(B-7) 


Consideration of possible nozzle exit geometries indicated that an exit- 
to-throat area ratio Ae/A=!' = 3. 5 and ratio dA/Aj = 1. 56 appeared practical 
within the constraints of the external geometry. This fixed M-/M =1.4 and 
T°/T° = 2.62 (nitrogen). Under these conditions, the ratio or stagnation 
pressure in the jet to that in the free stream is p9/p° = 3.468. 


Thus, in horizontal planes through the nozzle portion of the thermal fin, 
the following nozzle exit geometry was chosen: 


w 


w 

w e 

w w=i^ 
e 


S w>:< = 1. 56 (3. 5 )w- 
= 5. 46 w* 



d L ^ 

W W 

e 


The' thermal-fin Mach-plane area distribution in the region embodying'the 
nozzle is then given by 




^(y)dy = Kv^/^ 

s in/t ^ 
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where 


A = contrlbutLon of sol'id-fln area to Mach-plane area distribution 
w(y) = width of streamtube area due-to-heat as given above , 
y = distance in vertical direction 

' y. = value of y at initial point of heated-gas flow 
fJ. - Mach angle 

- 1/2 

K = constant = 0.00611 cm for thermal fin model 1.5 scale of 
i solid-fin model 

For w = 5.46 W''S U - 30 deg, hence 

+ 10.92w*(y - y.) = Ky^^^ (B-8) 
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APPENDIX C 

HEATING- TUBE ANALYSIS - 
DETERMINATION OF NOZZLE MASS -FLOW RATE 


From Appendix B., the design condition on the heater is that it delivers 
gas at the nozzle exit having stagnation temperature 

T° = 2.62 T° (C-1) 

At tunnel operating conditions, T° = 570°R (317°K). Therefore 

/ 

T° = I492°R (829°K) 

is the required jet stagnation temperature. 

The heating -tube design is based on the following analysis. Consider a 
tube of cross-sectional area S, diameter D, and wall temperature T with a 
gas at temperature T flowing through it. 



For heat transfer between the tube wall and the gas the convective heat- 
transfer coefficient is defined by 


q = hA(T - T, ) (C-2) 

^c w b 


where q^ is the rate of convective heat transfer to the gas in-the tube, h is the 
convective heat -transfer coefficient defined by the above relationship, A is the 
surface area of the tube, T^ is the wall temperature and T^ the bulk tempera- 
ture defined by 



r 

w 

pVc T ZTTrdr 
P 
o 


r 

w 

PVc 2m:dr 
P 
o 




(C-3) 
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where p is the gas density, V the velocity, c the specific heat at constant 
pressure, and r the radial variable. ^ 


For the case of constant -density flow in the tube (low-speed flow) and 
turbulent pipe flow for which 



(C-4) 


where subscript c denotes centerline values, there is obtained 



(C-5) 


Thus, to good approximation, ^ T^. 

To obtain the relationship between the gas temperature, the wall tem- 
perature, the power applied to the tube, and the tube electrical properties, 
consider the element 


t 



dx 


Then 

2 

q = h7TDdx(T - T) = ^ [1 +a(T - T„ )] dx -€aT^7TDdx (C-6) 

where I is the current, R the resistivity, a a coefficient giving the change of 
resistivity with temperature, Tp the reference temperature for resistivity, 

€ the wall emissivity, a the Stefan -Boltzmann constant, D and t the tube dia- 
meter and thickness, respectively, and the subscript c has been dropped on 
the temperature T without loss of generality. The first term on the extreme 
right-hand side of equation (C-6) represents the heating rate due to resistance 
and the second the rate of heat loss due to radiation. 

-4 

For the nichrome tubing under consideration, R = 10 ohm-cm 
(0. 3035 X 10-4 ohm. in. ) at 293°K (52?OR), Ot = 0. 0004/°K (0. 00022/OR), 

D = 0. 3175 cm (0. 125 in. ), t = 0. 0165 cm (0. 0065 in. ). At typical operating 
conditions, I = 60 amps, T — 1000°K, and the radiation energy is about 15 
percent of the total energy. However, the heating tube will be encased in a 



radiation shield and a portion of the radiant energy will be radiated back to it. 
Thus, neglecting the radiation term there is obtained 


T (TTDh -aP) - P(1 -aT„) 
w j R 

rrDh 


(C-7) 


2 

where P = I R/7rDt. 

Now, to obtain the gas temperature as a function of tube length, we 
make use of the fact that the wall temperature, gas temperature, tube mass- 
flow rate, and specific heat of the gas are related by 


me dT = 7rDh(T - T)dx 
p w 


(C-8) 


Substituting for T from equation (C-7) and integrating. 


there is obtained 


X 

D “ 



1 - 


ap 

77-Dh 


1 +Of(T - 


dT 


, (C-9) 


where Tq is the temperature at x = 0. The heat -transfer coefficient is given 
by (ref. 19) 


Nu 


D 


~= 0. 023 Re°* ®Pr°‘^ 
k D 


(C-10) 


where = PVD/n and Pr = c /i/k, with p being the gas density, V the velo- 

city, /i the viscosity, c^ the specific heat, and k the thermal conductivity. But 

PV = ^ (C-11) 


where S is the cross-sectional area of the tube. Taking the Prandtl number 
as 0. 7, there is obtained 


TTDh cs! 0, Ik 


0 . 2 


me 


D 


0 . 8 


(C-12) 


/ 

In addition, in the range 290°K < T < 900°K (522°R < T < 1620°R), k is given 
by 
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where P - 0.56, Icq = 2. 54 x 10~"^ (W /cm°K) [4. 07 x 10 ^ (B/ft-sec-°R )] . 
Defining M =mCp/PD and taking there is obtained from equations 

(G-9), (C-12) and (C-13) 


X 

D 


^ In [1 +a(T - T^)] - loot 


'MP 


o 


■) 




T-T 


dt 


(1 +ctt) 


(l, + ^t) 


0.2 


(C-14) 


where t s T - T^. Equation (C-14) gives the relationship between tube length 
and gas temperature, T, at the end of that length. Since the flow at that 
point will be subsonic, T — T°. 


Now, it is desirable to maximize the mass-flow rate in the tube at the 
design temperature to maximize the amount of fin represented by heat. From 
the definitions of M and P, it is seen that M oc m/I^. Hence it is desirable to 
run the heater at its highest current capability. Moreover, a heating-tube 
length of 152. 4 cm (60 in. ) is fixed by the length required through the sting 
and the amount through the model to the nozzle. For this length, the maxi- 
mum current is 60 amps. Then 


, . 7TC tm 

M _ . p 

“ I^Ra 


= 96. 18 X ro"^ cm-sec/kg (17. 14 x 10^ in. -sec/lbm) 


for nitrogen with c = 0. 25' cal/gm°K (0. 25 Btu/lbm°R) and a heater tube 
having properties ^given above. Thus, for x/D = 152.4 cm/0.3175 cm and 


T-T, 


538°K (970°R) 


m = 6. 52 gms/sec (0. 0144 Ibm/sec) 
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APPENDIX D 

DESCRIPTION OF CONTROL CONSOLE 


The controls for operation of the sonic boom test models are located in 
a self-contained unit shown in figure 30. A schematic diagram and the oper- 
ating sequence is shown on figures 31 and 32, What follows is a description 
of the various control functions involved in this unit. 


D. 1 Heater Power 


The power for the nichrome heater in the model is derived from 240 
volt, single phase, 60 cycle service. (208 volt service may be used with a 
corresponding derating of the unit, ) The line voltage is varied with an auto- 
transformer that controls the primary voltage on a two and one -half to one 
centertapped step-down transformer. This transformer feeds a full wave, 
silicon rectifier that is filtered with an incandescent light bulb and electro- 
lytic ‘capacitor. Output of 0-50 V dc and 0-100 amps is available at the re- ' 
ceptacle on the rear of the cabinet. An interlock is provided on the power 
system to ensure that the autotransformer is set to zero before the dc power 
can.be turned on. 


D. 2 Control Power 


Power for all relays and instrumentation is derived from the 240 V ac 
source through a step-down transformer and a fused 120 V ac circuit. No 
functions are operational until this circuit is activated as indicated by the 
white pitot light on the front panel. 


D. 3 Actuator Power 


Direct current for the motor which drives the model sting is obtained 
from a full wave bridge rectifier and a dropping rheostat (for speed control). 
Local setting of speed and direction is possible with controls on the front 
panel. Remote actuation (at full speed) is provided by relays that are ener- 
gized through a three -pin receptacle at the rear of the unit. (115 V ac relays 
are installed, but other voltages may be used by replacing K2 and K3 with 
octal base relays of appropriate coil voltages. ) 
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Figure 30. Sonic -Boom Control Console 
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D. 4 Gas Supply- 


Nitrogen for the heater is controlled by a throttling valve that operates 
on a gas supplied at a regulated, constant pressure. A sonic orifice is used 
for metering the flo-w so a pressure ratio of at least two must.be maintained 
at all times between "supply" and "heater” pressure. 


D. 5 Instrumentation 


All system parameters are measured and displayed on meters on the 
front panel. No recorder connections are provided. Power for transducers 
is furnished from a ±20 V dc, filtered, one-half wave rectifier. Pressure 
and position transducers are potentiometer devices. Nozzle temperature is 
measured with a chromel-alumel thermocouple, the output of which is ampli- 
fied in an integrated circuit unit for direct reading on the panel meter. 
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